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We propose physical interpretations for stochastic methods which have been developed recently to 
describe the evolution of a quantum system interacting with a reservoir. As opposed to the usual 
reduced density operator approach, which refers to ensemble averages, these methods deal with 
the dynamics of single realizations, and involve the solution of stochastic Schrodinger equations. 
These procedures have been shown to be completely equivalent to the master equation approach 
when ensemble averages are taken over many realizations. We show that these techniques are not 
only convenient mathematical tools for dissipative systems, but may actually correspond to concrete 
physical processes, for any temperature of the reservoir. We consider a mode of the electromagnetic 
field in a cavity interacting with a beam of two- or three-level atoms, the field mode playing the 
role of a small system and the atomic beam standing for a reservoir at finite temperature, the 
interaction between them being given by the Jaynes-Cummings model. We show that the evolution 
of the field states, under continuous monitoring of the state of the atoms which leave the cavity, 
can be described in terms of either the Monte Carlo Wave-Function (quantum jump) method or 
a stochastic Schrodinger equation, depending on the system configuration. We also show that the 
Monte Carlo Wave-Function approach leads, for finite temperatures, to localization into jumping 
Fock states, while the diffusion equation method leads to localization into states with a diffusing 
average photon number, which for sufficiently small temperatures are close approximations to mildly 
squeezed states. 

PACS numbers; 42.50.Lc, 42.50.Ar, 42.50.-p 



I. INTRODUCTION 

The dynamics of dissipative quantum systems is fre- 
quently described through a master equation for the re- 
duced density matrix, obtained by tracing out the degrees 
of freedom of the reservoir and making the Born-Markov 
approximation jl]. As usual in quantum mechanics, the 
corresponding solutions refer to ensembles of identical 
systems. In recent years, the attainment of low tem- 
peratures and low-dissipation regimes, as well as the im- 
provement of detection techniques, has allowed the inves- 
tigation of the dynamics of continuously monitored sin- 
gle quantum systems. Remarkable examples of these are 
single ions or Bose-Einstein condensates in electro- 
magnetic traps, probed by laser beams, and electromag- 
netic fields in high-Q cavities, probed by beams of highly 
excited atoms (Rydberg atoms) Q . This new generation 
of experiments, combined with the difficulties usually en- 
countered in solving the master equation, have stimu- 
lated the development of new techniques, which seek to 



describe the dynamics of quantum dissipative systems by 
stochastic evolutions of the corresponding state vectors 

IHl; 

Basically two approaches have been proposed, one 
which involves random finite discontinuities in the dy- 
namics of the system (Monte Carlo Wave-Function or 
"quantum jumps" method) |^,^ and another for which 
the stochastic evolution is generated by a diffusive term 
in the Schrodinger equation for the state vector 
These alternative techniques can be shown to be formally 
equivalent to the master equation approach , and in 
some cases do lead to a dynamical behavior resembling 
the experimental monitoring of a single realization. How- 
ever, they are often regarded as mathematical tricks, with 
no relation to a concrete physical evolution of the system. 
The advantage of using them, from the numerical point 
of view, is that one deals with state vectors, instead of 
density matrices, thus reducing the total amount of ma- 
trix elements to be calculated. In addition, they may 
provide insights into the behavior of dissipative systems. 
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In fact, and because of these two points, they have been 
extensively appHed to dissipative quantum systems, es- 
pecially in the fields of quantum optics ||,^^-|l|,^6[^] 
and solid state physics In some cases, these meth- 
ods have led to analytical descriptions of the dissipation 
process [^6 jl8|25]. More recently, the Monte Carlo Wave- 



Function (MCWF) method has been extended to non- 
Markovian interactions and to strong reservoir couplings 
beyond the Born and rotating- wave approximations pl| , 
and also to nonlinear master equations 13 1^. 

In this work, we show that it is possible to inter- 
pret the above-mentioned stochastic evolutions in terms 
of continuous measurements made on concrete physical 
systems, for any temperature of the reservoir. This pa- 
per extends therefore the pioneer work of Wiseman and 
Milburn [jll],^, who developed physical interpretations 
for Schrodinger stochastic equations describing the evolu- 
tion of a cavity mode in contact with a zero-temperature 
reservoir. They showed that the mode of the electro- 
magnetic field is described by a quantum jump equation 
if the outgoing light is directly detected by a photode- 
tector, while homodyne or heterodyne detection schemes 
lead to two different stochastic Schrodinger equations for 
the state of the field. 

We also consider the time-dependent behavior of a 
mode of the electromagnetic field in a cavity interact- 
ing with a reservoir (which may be associated to the 
continuum of modes of the field outside of the cavity, 
to which the internal field may be coupled via a semi- 
transparent mirror) . Our interpretations rely on the fact 
that the same master equation is obtained for quite dif- 
ferent reservoir models. We take thus as the reservoir an 
atomic beam which crosses the cavity, interacting reso- 
nantly with the field inside it. We show that the Monte 
Carlo or quantum jump method can be reproduced by 
taking as reservoir a beam of continuously monitored 
two-level atoms, prepared initially in a mixture of the 
two resonant states, in such a way that the population 
ratio between the excited and the ground state equals 
the Boltzman factor (see Fig. |^). On the other hand, 
a Schrodinger diffusion equation is obtained when the 
reservoir is assumed to be made of three-level atoms with 
a twofold-degenerated ground state, prepared in the same 
kind of statistical mixture as before, but now with the 
two ground state levels placed in a coherent superposi- 
tion. While a resonant exchange of energy is allowed 
between the cavity mode, the excited state and one of 
the ground state levels, the transition between the other 
ground state level and the excited state is assumed to be 
driven by an external (essentially classical) field (see Fig. 
1). 

In order to further explore the physics underlying these 
approaches, we calculate the evolution of an initial field 
in the cavity under continuous monitoring, for different 
types of initial states of the field, with the two different 
methods mentioned above. We find that for tempera- 
tures different from zero the two approaches lead to two 
kinds of localization in state space. For the quantum 



jump method, the state of the system approaches a Fock 
state, which suffers quantum jumps in such a way that 
the average distribution in time of the number of photons 
satisfies the thermal distribution. On the other hand, the 
diffusion equation leads to states which, for sufficiently 
low temperatures (average number of thermal photons 
smaller than one), are quite close to mildly squeezed 
states. Even though localized in phase space, these states 
have a diffusive behavior, spanning eventually a region 
of the phase space in such a way that again the time- 
averaged photon statistics coincides with the thermal dis- 
tribution. We will show that our atomic model for the 
reservoir allows a simple interpretation of these localiza- 
tion phenomena, which extends to finite temperatures 
the discussion made by Garraway and Knight [isl , and 
illustrates the general localization properties of quantum 
state diffusion equations demonstrated by Gisin and Per- 
cival 

In Section II, we review the stochastic approach to 
dissipative systems. In Section III we propose a phys- 
ical interpretation for the Monte Carlo quantum jump 
approach to the problem of field dissipation in cavity 
QED, for any reservoir temperature, while in Section IV 
we show how to interpret physically a description of the 
same problem based on a stochastic Schrodinger equa- 
tion. In Section V, we display our numerical results, and 
show that, depending on the physical procedure used to 
monitor continuously the field in the cavity, one may get 
localization in state space. Our conclusions are summa- 
rized in Section VI, while details of the calculations are 
displayed in the Appendices. 



II. STOCHASTIC SCHRODINGER EQUATIONS 
AND DISSIPATIVE SYSTEMS 

A wide class of master equations describing the evo- 
lution of dissipative quantum systems can be written in 
the Lindblad form [^: 



where 



C = Co + Cn , 

n 

2 

CoPs = t:[ps,Hs] , 



(1) 

(2) 
(3) 



^nPS = -^[ClCnPS + PsClCn] + CnPsC^ , (4) 

ps is the reduced density operator for the "small" system 
S (obtained by tracing out the degrees of freedom of the 
reservoir R from the density operator for the full system 
S+R), and Hs describes the Hamiltonian evolution of the 
small system S in the interaction picture. The operators 
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Cn act on the space of states of the small system S, and 
express the interaction of S with the reservoir R. The 
number of them depends on the nature of the problem. 

An example of such an equation is the master equation 
for a field in a lossy cavity, at temperature T, given in 
the interaction picture by 

+ r(l + n){apfa^ - ^a^ap/ - ^P/o^a) , (5) 

where a and a' are the photon annihilation and cre- 
ation operators, respectively, n is the average number 
of thermal photons, given by Planck's distribution, and 
r = 1/tcav, where icav is the damping time. In this case, 
one could set 



Ci = ^/T{^+n)a, C2 = Vf^a} . (6) 
A formal solution of Eq. (|^) is 

p(i) = exp(£t)p(0). (7) 

Let us define 



JnP — dnpC\ , 



(8) 



and write 



p{t) = exp 1^0^ + 51 + (Cn -Jn)]tj p(0) . (9) 



Note that 



(£„ - J„) PS = - i {ClCnPs + PsCiCn) . (10) 



Applying Dyson's expansion to Eq. (g), we get: 

= ^ / dim / dtm-1--- I dti 
m=0-^0 Jo 

{S(t ^m) Jn)S{tfji tm — l) 

n 

x...(^j„)5(ii)M0), (11) 



Each term in the above double sum can be considered 
as a quantum trajectory, the reduced density operator at 
time t being given by the sum over all possible quantum 
trajectories For each of these trajectories, Eq. ([T^ ) 
shows that the evolution of the system can be consid- 
ered as a succession of quantum jumps, associated to the 
operators J„, interspersed by smooth time evolutions, as- 
sociated with the operators S{t). The probability of each 
trajectory is given by the trace of the corresponding term 
in Eq. (PI). 

From Eqs. dTo) and (I2h, we can write: 



S{t)p = N{t)pN{t)^ , 



where 



N{t) = cxp 



(14) 



(15) 



Therefore, if p is a pure state, then S{t)p is also a pure 
state. The same is true for JnP, with J„ defined by Eq. 
(^). This implies that a pure state remains pure, when a 
single quantum trajectory is considered. Note also that 
the evolution between jumps is given by the non-unitary 
operator N{t). 

It is clear from Eq. (||) that different choices of the 
jump operators are possible. These different choices cor- 
respond to different decompositions in terms of quantum 
trajectories of the time evolution of the density opera- 
tor ps and, eventually, to different experimental schemes 
leading to the continuous monitoring of the evolution of 
the system. It is precisely due to this continuous mon- 
itoring that an initial pure state remains pure, since no 
information is lost in this situation: for a field in a cav- 
ity, this continuous monitoring amounts to accounting 
for every photon gained or lost by the field, due to its 
interaction with the reservoir. 

We will discuss now two different realizations of the 
reservoir, for a field in a cavity, which will lead to a Monte 
Carlo quantum jump equation, for the first realization, 
and to a Schrodinger equation with stochastic terms, for 
the second one. 



III. SIMULATION OF A MONTE CARLO SSE. 



where 



S{t) exp ■ 



Equation (0) may be rewritten in the following way: 



I't2 

dtm-1 ■ ■■ dti 
Jo 



m=0{ni} 
{S{t — tm)Jn,-nS{tm — ^m-l) 
X ... J„,5(tl)}p(0). 



We exhibit in this Section a physical realization of the 
Monte Carlo method. The corresponding experimental 

(12) scheme is shown in Fig. 0. A monokinetic atomic beam 
plays the role of a reservoir R and crosses a lossless cavity, 
interacting with one mode of the electromagnetic field. 
The cavity mode plays the role of a small system S. The 
atoms, regularly spaced along the atomic beam, are pre- 
pared in one of two Rydberg states: an upper state \a) 
or a lower state |6). The transition frequency uj between 
these two states is assumed to be resonant with the cav- 
ity mode. A similar model of reservoir was adopted in 

(13) Section 16.1 of Ref. ||. 
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The state of the atoms is measured by a detector just 
at the exit of the cavity. The ratio between the flux of 
upper state atoms ra and the lower state atoms Vb before 
their entrance into the cavity is chosen so that 
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(16) 



where huj is the difference in energy between \a) and 
and, as will be shown in the next paragraphs, T is the 
reservoir temperature. The constant ks represents the 
Boltzmann constant and n, given by Planck's formula 
(n = [exp{hw/kBT) — 1]^^), is the mean occupation of 
the modes with energy Hw in a bath at temperature T. 

We analyze now the time evolution of the state vector 
|\E'(t)) of S, under the continuous measurement of the 
atoms after they leave the cavity. We also assume that 
one knows the state of each atom before it interacts with 
the cavity. This may be achieved by selectively exciting 
the atoms to \a) or according to the proportion given 
by Eq. (p6|). We will adopt the following simplifying as- 
sumptions: (a) the atom-field interaction time r is the 
same for all atoms; (b) the spatial profile of the electric 
field is constant; (c) the cavity is perfect, i.e., the field 
state is changed only by the atoms; (d) the atom-field 
coupling constant A and the interaction time r are both 
small, so that the atomic state rotation is very small; 
(e) the rotating-wave and dipole approximations will be 
used; and (f), according to the statements (d) and (e), 
quantum cooperative effects will be neglected. In this 
case the interaction Hamiltonian in the interaction pic- 
ture will be: 



H = hX- 



(a a' 



\a){b\a) 



(17) 



The operators a and are annihilation and creation op- 
erators, acting on the space of states of the field mode. 
Just before the ith atom enters the cavity, the state de- 
scribing the combined system (atom i + field) is given 

by 



|M'a-/(t.)> = |*(t,)>®l*a(t.)). 



(18) 



Here |^a(ii)) — 1") \^a{ti)) — |&), depending on the 
state to which the atom was excited just prior to entering 
the cavity. 

At time U + t, the atom-field state vector, up to second 
order in r, is given by: 

\^a~f iU + t)) = (l- iXT\b) {a\ at - iXT\a) {b\ a 

- ^\b){b\a^a-^\a)(a\aaA \^a-f{U)), (19) 



where the tilde indicates that the state vector is not nor- 
malized. The expansion (|l|) should be very good in view 
of condition (d). We assume that {ra + rb)T < 1, so that 
there is at most one atom inside the cavity at each instant 
of time. After this atom exits the cavity and is detected. 



one of the following four cases will be realized: 

i. The atom enters the cavity in state \b) and is detected 
in the same state. In this case, according to Eq. ([is|), 
the state of S at time t = ti + t will be given by 

mU+r)) ^ (^1- ^a^a^^iU)). (20) 

ii. The atom enters the cavity in state \a) and it is de- 
tected in the same state \a). In this case, 

\^iU+r))=(^l-^aa^y^{U)). (21) 

in. The atom enters the cavity in the state |6) and it is 
detected in the state la). In this case, 



\^{U + t)) -iATa|*(i,))- 



(22) 



iv. The atom enters the cavity in the state |a) and it is 
detected in the state |&). Then, 



(23) 



Note that in the cases i and ii a small change in the 
state of "S" takes place, while in the cases Hi and iv a 
big change may happen (quantum jump). However, these 
last two cases are very rare, due to the small change of 
the atomic state during the interaction time. 

We consider now the change of j^*) from t to t + dt, 
where the time interval St is large enough so that many 
atoms go through the cavity during this time interval 
{ua ~ raSt ^ 1, Ub — rbSt ^ 1), and also much smaller 
than tcnv /n{n) , where (n) is the average number of pho- 
tons in the state. This last condition, as it will be seen 
later, implies that the probability of a quantum jump 
during St is very small. In most of the time intervals St 
the atoms will be detected at the same state they came 
in, since the transition probability is very small. The 
evolution of l^f) during these intervals will be given by: 



1 



2 



ata) |*(t)) 

HbX^T^ 



ata |*(t)). (24) 



This result does not depend on the ordering of the upper- 
state and lower-state atoms. We also note that in the 
interaction picture the state vector does not evolve when 
there is no atom inside the cavity, since the only source of 
field dissipation is the interaction with the atomic beam. 

Equation (|2^) displays the interesting property that 
the wave-function of the system (and, consequently, the 
mean energy) may change even when there is no exchange 
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of energy between the system and the measurement appa- 
ratus (represented by the atoms in the present case) . An 
easy way to understand this effect physically is to imagine 
that all atoms are sent into the cavity in the lower state, 
and are detected in the same state after exiting the cavity, 
for a given realization of the system, which starts with a 
coherent state in the cavity. Then, even though there is 
no exchange of energy between the atoms and the field in 
the cavity, as time evolves the ground state component of 
the initial state should also increase, since the results of 
the measurements lead to an increasing probability that 
there is a vacuum state in the cavity. In other words, 
the fact that there is no quantum jump, for that specific 
trajectory, provides us with information about the quan- 
tum state of the system, and this information leads to 
an evolution of the state. This is closely related to the 
quantum theory of continuous measurement |^,^ and 
also to quantum non-demolition measurement schemes 
proposed recently [p7| . This problem is also very similar 
to that of a Heisenberg microscope in which even the un- 
successful events of light scattering produce a change in 
the quantum-mechanical state of the particle . 
We introduce now the following definitions: 



2^2 



{n - ?'a)A r 



2^2 



1 + n 



■X't 



C2 



'Tna) . 



Using these definitions and relation dig), Eq. 
be rewritten in the following way; 



St 



1 ^ C'^C'„ 



(25) 

(26) 
dH) may 

(27) 



If an atom enters the cavity in state \a) and is detected 
in the state |6), the state vector of S suffers a "quantum 
jump" , and one photon is added to that system. On 
the other hand, a deexcitation in S occurs if an atom 
which entered in \b) is detected in the state \a). The 
probability of this events to occur may be calculated by 
using Eqs. ( p6| ) and ( p2| ) or (|2^); thus, the probability of 
an excitation (action of ) to occur between t and t + St 
is given by: 



Spi^St('^it)\ciC\\-^it)). 



(28) 



The probability of a deexcitation (action of a) during this 
time interval is: 



Sp2^St('^{t)\ClC2\'^it)). 



(29) 



The probabilities Spi and Sp2 are very low, so that the 
joint probability of having one excitation and one deexci- 
tation during the same time interval St is negligible. One 
may therefore write: 



I'^it + st)) = cr^c. 



'5Ni 



1 2" 5-/ 



1^(0) 



(30) 



where SNi and SN2 are equal to one or zero, with prob- 
abilities Spi and Sp2 for SNi and SN2 to be equal to 
one, respectively. This may be represented by writ- 
ing the statistical mean M{SNm) — {Cl^C,n)St. Also, 
SN„iSNn — SN„iSnm- One should note that the instants 
of time in which the quantum jumps occur during the 
time interval St are irrelevant, since the jump opera- 
tors can be commuted through the no-jump evolution, 
the commutation producing an overall phase which goes 
away upon renormalization of the state. This can be 
easily seen by rewriting the no-jump evolution, during a 
time interval Stj < St, as an exponential: 



St, 



St, 



exp 



0[(%) 



and using that 



(31) 



(32) 



where Ai = -{Sj /2)r{l + ri) and A2 = iStj/2)Tn. 

The results of the measurement may be simulated by 
picking random numbers. The state vector given by Eq. 
(pO|) may be normalized as follows: 



ib{t + St)) 



-.SNi 



Co 



-.SN2 



\fclch \fclC2 
(1 - 5N,){1 - SN2) (1 - Y E 

[l-StY,{ClCm) 



m)- (33) 



In the above equation, the first two terms represent the 
possible jumps, each normalized, as in the Monte Carlo 
method, and the last term is the no-jump evolution con- 
tribution, normalized with the corresponding prefactor 
that rules out the jumps. From Eq. (|3^ ) one gets for 
\dm) =\i^it + St))-\yjit))- 



dm) = { E 



Cin Cry 



- 1 



SN„ 



I 5](Ct„C,„ - {ClCrn))\ I m)- (34) 

m ) 



IV. SIMULATION OF THE HOMODYNE SSDE. 

We show now that, by a suitable modification of the 
atomic configuration, it is also possible to interpret phys- 
ically the stochastic Schrodinger equations in terms of 
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continuous measurements made on atoms which cross the 
cavity containing the field. The corresponding scheme is 
shown in Fig. 2: a beam of three- level atoms with a de- 
generate lower state (states b and c) crosses the cavity, 
the field in the cavity being resonant with a transition 
between one of the two lower levels (say, level b) and the 
upper atomic state a, while a strong classical field con- 
nects the other lower state with the upper level (one may 
assume that both fields are circularly polarized, so that 
the cavity field cannot connect a and c, while the strong 
field does not induce transitions between a and b). 

We also assume that the atom is prepared in either a 
coherent superposition of the two lower levels: 



i'aton.) = ^(1 b)+ I C)), 



(35) 



or in the upper one, following a Boltzmann distribution 
corresponding to a temperature T for the atoms, which 
act as a reservoir for the quantum field in the cavity. 
In the interaction picture, one can write: 



H = hgacie I a)(c | +s \ c){a |) 
+ hgab{a^ \b){a\+a\ a){b \ 



(36) 



We assume for simplicity that gac — gab = <?, and that 
e is real. The time evolution of the wave function , to 
second order in the coupling constant is: 



1 - 



iHt H'r' 



2h^ 



\m)- (37) 



As in the previous model, there are two possible quan- 
tum jump processes. The first one corresponds to the 
atom entering the cavity in the coherent superposition 
of lower states, and being detected in the upper state. 
After the measurement, the state of the field is given by: 



{e + a)\Mt))f- 



(38) 



The corresponding probability of detecting an atom 
in I a), after a time interval St, staring from the initial 
superposition state, is given by: 

2 2 

Spi = nb^i^bfit) I (£ + at)(e + a) \ ^Pf{t)), (39) 

where nb = rbSt, rb being the rate of atoms injected in 
the superposition of the lower states. 

The second jump process corresponds to the atom en- 
tering the cavity in the upper state |a), and being de- 
tected in the superposition of lower states. Then, the 
state of the field after the measurement is: 



ie + a^)\m)f- (40) 



The corresponding probability is given by: 



2 2 

SP2 = na^ii^fit) \ie + a){e + a^)\ i^f{t)) , (41) 

where Ua = radt is the number of atoms which enter the 
cavity in state |a), during the time interval St. 

This analysis suggests that the quantum jump opera- 
tors corresponding to these two processes should be, re- 
spectively. 



and 



where 



Ci = Vr(l + n)(e + a) 
6*2 = Vf^{e + flt), 



in - To) 



g^r^ 



2 2 

n a t"- 
\ + n 2 



2 2 

n 2 ' 



(42a) 
(42b) 

(43) 



Formally, these jump operators arc retrieved by rewrit- 
ing the Master Equation (^) in the following equivalent 
form: 

^ = ( Ji + J2)P} - ^-^^ [(ata + 2ea + e^)pf 
Tn 

+ 2ea + e^)] , (44) 

with: 



i = l,2 



(45) 



being associated with the jumps, the operators Ci being 
now given by Eqs. (^2[). 

We derive now the stochastic Schrodinger equation 
that describes the present measurement scheme. 

With the above jump operators, and using the expan- 
sion given by Eq. (|l^), we show in the appendix A that 
the joint probability of getting mi and TO2 jumps corre- 
sponding respectively to the first and second processes 
described above is given by the following expression: 



.(At) 



exp/xi- 

mi 

X Trjexp^' 
1 



exp fi2 - 



{P2) 



TO2! 



1 H — (mifl + 17120^) 

e 



1 H — (mia^ -t- m2a) 



exp (3^ 



(46) 



where: 



Ail =rAte2(l + n), (47) 
Ai2 =rAte2(7i), 
TAt 

13' = — [a)a{2ri + 2ea(n + I) + 2ea)n + n)] . 



6 



From Eqs. ( ^6| ) and (^), one can readily find {rrii) and 



Up to order £~^", one finds: 



mi 



Mi 



with 



(48) 



(49) 



Going back to the definition of S{t), one may write: 

S{M) = N{At)pN\At), (50) 

in terms of a smooth evolution operator N that preserves 
pure states. This operator N is given by EclJM). with 
the jump operators Cm now given by Eqs. (42|)7Now, if 
we consider a sequence of jumps (of the two kinds, in the 
present analysis) and evolutions, the state vector of the 
field will evolve according to: 



On the other hand, during the "no-click" periods, the 
nature of the coherent state is preserved, changing only 
the coherent amplitude, all the way to the vacuum. 

This situation has been previously studied |^ in the 
context of the continuous measurement theory of three- 
level atoms and two resonant fields, with the difference 
that in that work the number of detections was a prede- 
termined quantity. However, the net result of the preser- 
vation of the coherent nature of the state of the field, 
along the trajectory, is the same. 

b) If we start with a Fock state, the quantum jumps 
will invariably produce a mixture of various Fock states, 
while the waiting or "no-click" periods will only generate 
numerical factors in front of those Fock states. 

In the finite temperature case, the situation is more 
complex, since there will be also creation of photons, that 
will disturb an initial coherent state and produce further 
mixtures in the Fock state case. 

A more detailed analysis of these various cases is de- 
scribed in the next section, devoted to the numerical sim- 
ulation. 



I ij)f{At)=N{At~tm)C2N{U^ 
= N{At)C:^^C'^' I 7/.)/(0). 



<,„_i)Ci... I V)^(O) 
(51) 



In the last step, in deriving Eq. (pl|), we used that the 
commutators between the jump operators and the no- 
jump evolution produce overall phases, like in the Monte 
Carlo evolution given by Eq. (|30|). 

Now, we consider mi, i = 1, 2 as a couple of random 
variables with non-zero average, and write them as: 



V. NUMERICAL RESULTS AND 
LOCALIZATION. 

We present now the numerical calculations correspond- 
ing to the two equations associated with the two mea- 
surement schemes discussed above. We consider in these 
calculations the general case in which the temperature of 
the reservoir is taken as different from zero. 



{rrii) 



AW, 



(52) 



where the AWi are two real and independent Wiener 
increments, with: 



{AWf 



At, 



1,2. 



(53) 



From Eqs. (|T|) and (|5^ ) and up to order e we 
get the following Homodyne Stochastic Schrodinger Dif- 
ferential Equation (HSSDE): 



A" 



I i>)f{At) =\ ^)f{At)- I V)/(0) 

r r 

--{\+n)a^a - -{n)aa) + 2T{Xi){a{l + n) 



a)n)\ At - 



iv/r(l + n)AIFi| I i>)f{Q). 



(54) 



At zero temperature, a typical quantum trajectory in 
this homodyne scheme is as follows: 

a) If one starts from a coherent state, the quantum 
jumps will only produce a multiplicative factor in the 
wave function of the field, factor that can be absorbed in 
the normalization. 



A. Quantum jumps evolution 

We consider first an example in which the initial state 
of the system is a Fock state with three photons. We 
assume that the temperature of the reservoir corresponds 
to an average number of photons also equal to three. The 
corresponding evolutions is exhibited in Fig. 3. The state 
of the system remains a Fock state, with a number of 
photons which keeps jumping between several values, in 
such a way that the average number of photons is equal to 
three. We have verified that the probability distribution 
for the number of photons is a Bose-Einstein distribution, 
as long as the observation is done over a sufficiently large 
time. 

Figure 4 displays two different views of the evolution of 
the photon number population |a„p of an initial coher- 
ent state. These figures clearly exhibit the dual nature of 
the system dynamics, with quantum jumps interspersed 
by non-unitary evolutions. In the displayed realization, 
the vacuum component of the state increases until the 
first quantum jump occurs. This jump corresponds to 
the addition of a thermal photon to the system, lead- 
ing to the disappearance of the vacuum component. The 
second jump corresponds to the absorption of a photon 
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from the cavity field, leading to the reappearance of the 
vacuum state. The combination of the non- unitary evolu- 
tion with the quantum jumps finally leads to a Fock state, 
which under the action of the reservoir keeps jumping, in 
such a way that the photon number distribution over a 
long time span reproduces the Bose- Einstein distribution. 
This process is illustrated in Fig. 5, which displays the 
time evolution of the Q distribution for the field, defined 
for each realization as Q = \ {a\'4))\'^ / n , where \a) is a co- 
herent state with amplitude a. The initial Q distribution 
is a Gaussian, corresponding to the initial coherent state 
|ao); with ao — -\/15/2(l -f i). This distribution evolves 
into the one corresponding to a Fock state, with a num- 
ber of photons which keeps jumping around the thermal 
value rl = 2, in the same way as shown in Fig. 3. 



B. Diffusion-like evolution 

We consider now the evolution corresponding to the 
situation displayed in Fig. 2. We consider as initial state 
the same coherent state as in Fig. 5, the reservoir temper- 
ature being also the same as before {n — 2). In this case, 
the system evolves according to the homodyne stochastic 
Schrodinger equation given by Eq. ([54[). After some time, 
the Q function approaches a distorted Gaussian, with a 
mild amount of squeezing along the direction of the axis 
corresponding to the real part of a. The center of this 
Gaussian keeps diffusing in phase space, so that after a 
long time span the time-averaged distribution coincides 
with the Bosc-Einstcin distribution. Similar localization 
patterns were demonstrated in Refs. |3Q,pi[. 



C. Analytical proof of localization 

For the quantum jump situation, it is actually possi- 
ble to demonstrate that the system evolves in the mean 
towards a Fock state, for non-zero temperatures. 

We first define two kind of variances, for an arbitrary 
operator O. 

For the Hermitian case: 



and for the non-Hermitian case: 



(55) 



I AOp =(Ot-(Ot))(o_(o)) 

= Oto- (Ot)o_ot(o) - (Ot)(o), (56) 

so that 

(I AOp) = (Oto)_(ot)(o). (57) 
In particular, we are interested in two quantities: 



Qi = (l Aap), 
Q2 = (| Anp), 



(58) 
(59) 



that measure the distance of the state from being a co- 
herent or a Fock state, respectively. 

We start with the quantum jump equation: 



\dip) =--H\ ijj)dt 



\Y.{ClC,r.~{Cl){C„,))\^)dt 

m 

Cm. 



™ \l CmCm 



1) I ^)5N„ 



with: 



M{5N^) = {ClC,n)dt, 



(60) 



(61) 
(62) 



We will calculate, using Ito's rule of calculus, Qi 
a nd Qa fo r T = (C = yTa) and T > (Ci = 
V(n + l)ra, C2 = VT^at). 

We first develop some general expressions, which will 
be applied to calculate the above variances. 



d (O) = (d?M O I V) + (?M O I di;) + {dtfj \ O \ dij^) 
^-'-{[0,H])dt~\{{0, 
((CtOC) - (CtC)(0)) 



-^-{[0,H])dt ~]^{{0.,C^ C})dt + {0){C^ C)dt 



(CtC) 



-5N, 



(63) 



and similarly for the case in which several jump operators 
are present. 

For the variance of a non-Hermitian operator, we have: 

d{{\ AO H) = d{0^0) - {0)d{0^) - {0^)d{0) 

- d{0^)d{0). (64) 

After a simple calculation, one gets: 

d{{\AO = AO\',H])dt 

-^{{\AO\',C^C})dt 

+ (I AO f)(C^C)dt-{\ AO p)5iV 
{C''0''OC){C''C) - (C^O^C){C''OC) 



+ 



{C1C){CW) 



SN. 



(65) 



In the Hermitian case, on the other hand, we get: 



^{[AO\H])dt~^{{AO\C^C})dt 



d ((AO^)) . 

+ {AO^){C^C)dt~ (AO^)SN 

{C'<O^C){C''C) - {C^OC){C^OC) 



{CW){CW) 



SN. 



(66) 



Now we specialize to several cases: 

a,) T — 0, O = a, C = vTa, and H — Tiujo^a. 

Using the above general expressions, we write: 
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Aa p)) = [-r(ataata) _ 2r(a^a) (a^) (a) 

r r 

+ T{a^a){a'^a) + - {a'^ a'^ a) (a) + -{a^aa''){a) 

r r 

+ -{aa''a){a'i) + - {a'' aa) {a^)]dt 



{a''a)5N + {a''){a)5N 
{a^ a^ aa) {a^ a) — {a^ a^ a) {a'' aa) 
(a'I'a) (at a) 



SN. 



(67) 



The above results are neither strictly positive or nega- 
tive, so we cannot draw any conclusion; however, for the 
statistical mean: 



M'M^ = .ri\Aa\i^) 



T({Aa^)a^a)(a^aAa) 



<0, 



(68) 



so, in the mean, the system goes to a coherent state, 
which, in this case, is the vacuum. 

b) T > 0, O = a, Ci = ^(n + l)ra, C2 = VT^a'^ , and 
H = huja^a. 

The reader can easily verify, with a little algebra, that, 
in this case, neither d{{\ Aa p)) or Md{{\ Aa p)) are 
strictly negative. 

c) T > 0, O = atfl, Ci = v^(n + l)ra, C2 = VT^at, 
and 7? = fiLoa'^a. 

In this case, as shown in appendix B, d{{Aa'^ a)^) is not 
negative, but M d{{Aa'^ a)'^) is: 

^^ d{{Aa^af) ^ ((Aata)ata)(at„(Aat„)) 



dt 



r(n) 



(atfl) 

((Aaat)aat)(aat(Aaat)) 



<0. (69) 



So is strictly diminishing in the mean, even at 
T > 0. Since Qi is not, the final state will not neces- 
sarily be the vacuum. Indeed, there is no unique final 
state in this case. It is easy to show from Eq. (|6^) that 
M[d{{Aa'^ a)"^) / dt] = if and only if the state of the sys- 
tem is a Fock state. This result shows therefore that any 
initial state approaches eventually a Fock state | n) , with 
n fluctuating with a mean n. While the exact relation be- 
tween the ensemble average behavior and the long-time 
behavior of a single trajectory is not completely obvious, 
it is clear that the probability of a trajectory violating 
these inequalities over a long period is very small. Once 
the trajectory approaches a Fock state, it remains a Fock 
state for all time thereafter. This is reflected in the nu- 
merical results. 



the damping of one fleld mode in a cavity at temperature 
T. 

This field damping mechanism can be modeled as an 
atomic beam, whose upper and lower population ratio is 
given by the Boltzmann factor, crossing a lossless cavity. 

The quantum jump trajectory can be interpreted as a 
continuous monitoring of the outgoing two-level atoms, 
which are resonant with the cavity mode. We show both 
numerically and analytically that this continuous mea- 
surement on the reservoir leads, for each trajectory, to a 
pure Fock state. At a later time and due to the non zero 
temperature, a thermal photon may produce a jump to 
a different Fock state, thus leading, as time goes on, to a 
series of Fock states, whose statistics will reproduce the 
thermal distribution. 

In the case of the Homodyne Stochastic Schrodinger 
Differential Equation, the proposed damping mechanism 
consists of a three-level atomic beam, with a split ground 
state, whose population ratio of the upper and lower lev- 
els is given by the Boltzmann factor. The atoms cross 
again a lossless cavity, being resonant with the mode of 
the field under consideration. A second field is externally 
applied, with the same frequency but different polariza- 
tion, so that each of the two fields connects the upper 
atomic state with a different lower sub-level. If this ex- 
ternal field is a strong classical field, we show analytically 
that the stochastic Schrodinger equation describing the 
behavior of the quantum field in the cavity corresponds 
precisely to the Homodyne Stochastic Schrodinger Dif- 
ferential Equation. 

The beam is then continuously monitored as it exits 
the cavity. Numerically, one observes, for low tempera- 
tures, that the state of the field goes to a mildly squeezed 
state, centered around a value of a which diffuses in phase 
space, in such a way that the time-averaged distribution 
again reproduces the thermal state. 
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VI. CONCLUSIONS. 



APPENDIX A: DERIVATION OF THE HSSDE. 



We propose here a physical interpretation of the 
Quantum Jump approach and the Homodyne Stochastic 
Schrodinger Differential Equation, using as an example 



Here we present the detailed derivation of the Homo- 
dyne Stochastic Schrodinger Differential Equation. 
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We start from the expansion given by Eq. (|11|), which 
in the two-jump situation, and neglecting the commuta- 
tors between the jump operators and the no-jump evo- 
lution (for the same reason as discussed in the previous 
section), can be expressed as: 



(At) 



mi +rn2 



mi\m2\ 



-SiAt)JpJ^'p{0). (Al) 



The probability of mi and TO2 quantum jumps of the 
respective types, is given by: 



mi\m2\ 



(A2) 



The Master Equation of the field, corresponding to a 
lossy cavity at Temperature T, may be written as: 

^ = ( Ji + J2)p - [a^a{l + 2n) + 2e(l + n)a^ 



+ 2ena. + n + e'^{l + 2n)] 
- - [a^a{l + 2n) + 2e(l -I- n)a + 2ena^ 
+ n + e'^{l + 2n)] p. 



(A3) 



Therefore, according to the discussion in Section II, 
one possible way of writing S{ At) is: 



S{At)p:^ N{At)pN{At)\ 



(A4) 



with: 



N{At) = exp [a^a{l + 2n) + 2e(l + n)o) 

2ena-f n-f £2(l + 2rT)] |. (A5) 
Using Eqs. ([A^) and (|A5|), we can write: 



Pyni ,m2 ( ^0 

Tr 



exp^i(/ii)'"i 



mi 



exp P2{p2y 



TO2 



exp(/3')(l + -T'' (1 + -)"V(1 + —T' 

£ £ e 

a , 



X (1 + -)"^ exp(/3t') 

£ 

where: 

pi = VAt£^[l + n), 
P2 TAte^n, 
— FAi 

/?' = — {a)a{l + 2n) + 2 [e(l + n)a + £na^] + n] . 

From Eq. ([A^), we can now calculate (m,j) and (7^^ = 



(A6) 



(A7) 



TO, 



(wi) up to order (-) ^ • The result is: 



(m,) =^,(l + -(Xi), 
e 

erf = Alj. 



(A8) 



Now, we turn to the final step of this calculation, which 
yields the time evolution of the state vector. 

After repeated jumps and no-jump events, the unnor- 
malized wave function for the field can be written as: 

I 4>)f{At)=N{At~t^)C2N{t^-Un^i)CiN.. I V)/(0), 

or, except for an overall phase factor: 

I ^)f{At) = N{At)C^^Cr I V')/(0), (A9) 

where the tilde (^) indicates that the state vector is not 
normalized. 



Using Eqs. ( |A5| ) and (|A9D , one can write, up to a nor- 
malization constant: 

I V^)/(At)=exp L£M{ata(i+w) 



2 [£(1 + n)at + £na] } 

+ \ ™2 



X I 1 + — 



or, expanding, up to £ "^/^ 
TAt 



1 + ^)"' IV')/(0), (AlO) 



^>/(A<) = 



1 - 



{a) a{\ + n) + aa)n) 



FAte(a(l + n) + a^n) 

I, 

1 H — (toiC + m2a' ) 

£ 



V')/(0). (All) 



We are interested in the e — > cxd limit. In deriving 
Eq. ( All ) we considered £ large, TAt ^ e^'^/^, and 

mi,m2,Ati,M2 ^ e^^^- 

Now, we consider two random numbers with non-zero 
average toi and m2: 



TOi = (mi) -|- 
TO2 = (m2) + 



'At 

0-2 



AT^i, 
AW2, 



(A12) 



which satisfy: 

{{AWiY) = ((AT4^2)') = At. 



(A13) 



We notice that AWi are two independent Wiener pro- 
cesses. 



Finally, Eq. (All) can be written as: 
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I V)/(At) =1 V)/(Ai)- I V')/(0) 
-^{l+n)a^a - ^{n)aa^ + 2r(Xi)(a(l + n) 

+ a^/T{l + n)AWi | | ^)/(0). (A14) 
which is the desired result. 

APPENDIX B: FLUCTUATIONS. 

We want to calculate d{{Aa^a)'^) and Md{{Aa'^a)'^). 
We do it first in a simple case T = 0, O = a^a, C = 
VTa, and H = hua^a. 

d {{Aa^'af) = rSt{-{a'^aa'^aa'^a) + 2(a''aa^a){a^a) 

— 2(0^ a){a^ a){a^ a) + {a^ aa^ a) {a^ a)} 

- {a^aa''a)SN + {a''a){a''a)SN 
{a^ aa^ aa) {a^ a) — {a^ a"^ aa) {a^^ aa) 



+ 



{a^a){a'^a) 



-6N, (Bl) 



or: 



d ((Aata)2) = -T5t{{Aa^a){Aa'^a){Aa^a)) 
- {{Aa^af)6N 

{a^ aa^ aa) {a^ a) — {a^a^aa}{a^a'^aa} 
{a<a){a'a} 

Now, we apply the above resu lts to the more interesting 
case T > 0, O = aU, Ci = ^/{n + l)ra, C2 = VT^a^ , 

H = huja^a: 

d {{Aahif) = -T{n + l){{Aa^a)){Aa^a){Aa^a))dt 

- {{Aa^af)SNi 

^ {{a^ aa^ aa^ a) {a^ a) — {a'^aa)a){a^aa)a))5Ni 

{a t a){a^a) 
+ Tndt[-- {aa^ aa^ aa'\) + 2{aa^aa^) — (aa^) 
+ 2{aa^aa)){a^a) - 2{aa^){a^a) - {aa)){a) a){a^ a) 
+ {a^ ao) a) {aa^) — {a^ a) {a) a) {aa^)] 

- {{AaUf)8N2 
{{aa^aa^aa^){aa^) — {aa)aa)){aa^aa)))5N2 



+ 



{aa]){aa^) 



(B3) 



In the above expression, neither the deterministic or 
the stochastic term is definitely non-increasing. But in 
the mean it does decrease: 



M 



d{{AaUf) _ ^ ^ ((Aata)ata)(ata(Aata)) 



dt 



- rn 



[a] a) 

{{Aaa})aa^){aa\Aaa^)) 
(aaf) 



<0. (B4) 
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FIG. 1. Physical realization of a quantum jump trajectory. 
A beam of two-level atoms crosses a resonant cavity. 

FIG. 2. Physical realization of the homodyne stochastic 
Schodinger trajectory. A beam of three-level atoms crosses a 
resonant cavity, being subjected to an external classical field. 

FIG. 3. Quantum jumps for an initial Fock state with 
n = 3, the number of thermal photons being also equal to 
three. The setup is the one shown in Fig. 1 

FIG. 4. Two views (a and 6) of the evolution of an initial 
coherent state (average photon number equal to three) , in the 
quantum jump approach. The temperature of the reservoir 
corresponds to a number of thermal photons equal to 0.2. At 
Tt = 1.52 a photon is absorbed by the cavity mode, while 
around Pt = 3 a photon is lost by the field in the cavity. 
Before the first jump, the amplitude of the coherent state 
decreases exponentially. After some jumps, the state becomes 
a jumping Fock state. 

FIG. 5. Evolution of the Q function, for the quan- 
tum ju mp ap proach, and an initial coherent state, with 
ao = ■\/l5/2(l -|- i). The temperature of the reservoir cor- 
responds to a number of thermal photons equal to two. The 

initial Gaussian, corresponding to a coherent state, evolves 
into the distribution corresponding to a jumping Fock state. 

FIG. 6. Evolution of the Q function, for the d iffusiv e evo- 
lution, and an initial coherent state, with qo = ^ 15/2(1 + j). 
The temperature of the reservoir is the same as in Fig. 5. The 
initial Gaussian, corresponding to a coherent state, evolves 
into a distorted Gaussian, whose center diffuses in phase 
space. 
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